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Abstract 
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CN , In the leptophilic model, one Higgs doublet couples to quarks and another couples to leptons. 
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We study the supersymmetric version of this model, concentrating on the tightly constrained Higgs 
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sector, which has four doublets. Constraints from perturbativity, unitarity, and LEP bounds are 
considered. It is found that the lightest Higgs, h, can have a mass well below 114 GeV, and 
for masses below 100 GeV will have a substantially enhanced branching ratio into r pairs. For 
this region of parameter space, traditional production mechanisms (Higgs-strahlung, W fusion and 
gluon fusion) are suppressed, but it may be produced in the decay of heavier particles. The second 
lightest Higgs has a mass of approximately 110 GeV for virtually all of parameter space, with 



J> ' Standard Model couplings, and thus an increase of a few GeV in the current lower bound on the 

Standard Model Higgs mass would rule out the model. The two heavier Higgs are both gauge- 



phobic, one decays almost entirely into bb and can be produced via gluon fusion while the other 
decays almost entirely into r + r~ but can't be easily produced. 



* grmarshall@wm.edu 
i mtsher@wm.edu 



1 



I. INTRODUCTION 



The main purpose of the Large Hadron Collider (LHC) is the study of the mechanism of 
electroweak symmetry breaking (EWSB). One of the simplest and most studied extensions of 
the Standard Model is the Two Higgs Doublet Model (2HDM), in which two scalar doublets 
are jointly responsible for electroweak symmetry breaking and fermion mass acquisition jl, ^J. 
This model has a very rich phenomenology, including charged scalars and pseudoscalars. 
Among the earliest motivations for the 2HDM is its additional CP violation relative to 
the Standard Mode, ffl, which can provide aa additional mm of baryogenisis and the 



12| will 



relative abundance of matter to antimatter in the universe [10|, [ill]. It was also motivated 
by the fact that supersymmetric models and models with a Peccei-Quinn symmetry 
always require a minimum of two Higgs doublets. 

In order to avoid unobserved tree- level flavor changing neutral currents (FCNCs), all 



fermions with the same quantum numbers (and which are thus capable o: 



couple to the same Higgs multiplet. The Glashow- Weinberg theorem [l3[ states that a 



mixing) must 



necessary and sufficient condition for the absence of FNCNs at tree-level is that all fermions 
of a given charge and helicity transform according to the same irreducible representation of 
SU(2), correspond to the same eigenvalue of T3, and that a basis exists in which they receive 
their contributions in the mass matrix from a single source. In the 2HDM, this is due to the 
introduction of discrete or continuous symmetries. Generally one may either take both up 
and down type quarks to couple to the same doublet or have each couple to its own doublet. 
It is usually assumed that the leptons couple to the same doublet as the down type quarks, 
in which case the former scenario describes the Type I 2HDM while the latter describes the 
Type II 2HDM. Such couplings can be enforced by imposing a suitable Z2 symmetry, which 
may simply be imposed ad hoc or which may arise as a subgroup of a continuous symmetry 
(as in Peccei-Quinn or supersymmetric models). 

Despite the traditional convention that leptons couple to the same doublet as the down 
type quarks, there is no a priori reason why this must be the case. An alternative possibility 
is that both the up and down type quarks couple to one doublet while the leptons couple to 
the remaining doublet. While the traditional 2HDMs have received a great deal of attention, 
relatively little work has been done in investigating this alternative possibility. Those who 
have focused on this model 14j-|l8| have referred to it by several names, our selection of 
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which is the Leptophilic Two Higgs Doublet Model (L2HDM). As noted by Su and Thomas 
[jjj], the consequences of a L2HDM could drastically alter the possible detection channels for 
a light Higgs at the LHC, so it is important that it be considered as incoming data begins to 
arrive. Furthermore, the possibility of substantially enhanced leptonic couplings (which can 
only occur in leptophilic models) may shed some insight into explaining recent experimental 

16] . 



results from PAMELA, Fermi LAT, and H.E.S.S. 

There also remain alternative possibilities. One can couple the up-type quarks and leptons 
to one Hig gs doublet and the down-type quarks to the other (referred to as the "flipped" 



model 



191 ] ) or one can couple all of the charged fermions to one doublet and the right-handed 



neutrino to another (referred to as the "neutrino-specific" model) [20|. While interesting in 
their own right, these models do not offer the possibility of substantially enhanced leptonic 
couplings, and we will not focus on them. 

The most popular extension of the Standard Model is supersymmetry, which can solve the 
hierarchy problem and which has a very tightly constrained Higgs sector. Thus, one is led 
to consider the supersymmetric versions of these alternative 2HDM models. Recently, with 
McCaskey, we considered 21[ the supersymmetric version of the "neutrino-specific" model, 
and found some remarkable signatures, including pentalepton and hexalepton events with 
very high rates at the Tevatron and the LHC. In this work, we extend the L2HDM to in- 
corporate supersymmetry. The resulting Supersymmetric Leptophilic Higgs Model (SLHM) 
leads to exciting phenomenological prospects. In the scalar sector, the strong constraints on 
the Higgs potential will substantially alter the phenomenology of the lightest Higgs boson, 
since decays to leptons can be substantially enhanced, and the decrease in the coupling 
to the gauge bosons means that the current LEP bounds will not apply, and much lighter 
Higgs bosons can be tolerated. In addition, the supersymmetric partners to the leptons and 
the leptonic Higgs doublet are influenced by the unusual Yukawa structure. In the case 
of R-parity violation, the lightest supersymmetric particle (LSP) could decay into leptons. 
Without R-parity violation the LSP might annihilate into leptons [16]. In this paper, we 
will focus on the scalar sector, since the results may be testable in the very near future at 
the Tevatron. 

The layout of this paper is as follows. In Section |TT] we review the setup of the L2HDM. In 
Section [TTT1 we introduce the SLHM and calculate the scalar mass matrices. In Section [TVl we 
consider various constraints on the model's parameter space by focusing on the neutral scalar 
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sector. By combining results from Yukawa coupling perturbativity considerations, unitarity 
requirements, and direct searches for Higgs bosons at LEP, we obtain severe restrictions on 
the model's parameter space. In Section |V] we discuss the phenomenology of the lightest 
and next-to-lightest Higgs bosons at the Tevatron and the LHC, and then in Section IVIj we 
conclude. 



II. THE LEPTOPHILIC TWO HIGGS DOUBLET MODEL 



The L2HDM contains two scalar SU(2)l doubles Q q and <3?^. A discrete Z 2 symmetry is 
imposed under which — > — <Ev and e Ri — > — e Rt , but all other fields are invariant. The 
resulting Yukawa lagrangian is given by 



C Y = -{Y£u Ri & q ■ Q L] + Y*d Ri & q ■ Q L . + Y^\ ■ E Lj + h.c.} ; 



(1) 



where 



Qu 



and $x 



for X = q, £ and <J> g = io"2$ g . The Higgs sector potential is given by [14j, |22] 
V = m 2 q \$ q \ 2 + m £ 2 |$ £ | 2 + (m 2 ,$t$, + h.c.) + y |$ 9 | 4 + y \<Z> e 



+ A 3 |$ g | 2 |^| 2 



+ A 4 |$;$£i 2 + y[($;$£) 2 +h. c . 



(2) 



The physical scalars consist of two neutral scalars h and H, a pseudoscalar x°, and a 
charged pair H ± . The other three degrees of freedom are the Goldstone bosons G and G°, 
which are eaten by the W and Z° respectively. If one defines the mixing angle tan /3 = v q /ve, 
the physical charged scalars can be expressed as 




g 



(3) 



cos j3 sin j3 
— sin (3 cos (3 

The physical neutral scalar states are expressed in terms of the mixing angle tana, which 
can be solved for in terms of the entries of the neutral scalar mass-squared matrix tan 2a = 
2M 1 2 2 /(M 1 2 1 - M| 2 ). One then finds the following relation 




V2 



cos a sin a 
— sin a cos a 




(4) 



The vertex factors for the couplings between the charged scalar and fermions are given by 



15| 



tt+ ( ^tan/3 \ ^ 

For large tan /3 the neutrino-lepton coupling to H + is magnified while the quarks' coupling 
to H + is diminished. The neutral scalar couplings to the charged leptons will similarly be 
magnified. An interesting feature of the model is that tan can be much larger than in the 
conventional 2HDMs without causing problems with perturbativity and unitarity, since the 
Standard Model leptonic couplings are smaller than the quark couplings. 



III. THE SUPERSYMMETRIC LEPTOPHILIC HIGGS MODEL 

In this section we introduce the minimal leptophilic model required to incorporate su- 
persymmetry. A SLHM will require a minimum of four Higgs doublets in order to achieve 
anomaly cancelation. Therefore, we add to the MSSM two Higgs doublets H and Hi with 
weak hypercharge assignments +1/2 and —1/2 respectively. The four Higgs doublets along 
with their weak hypercharges are listed in the table. 







H d 
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H t 


Uy(l) 


+ 1/2 


-1/2 


+ 1/2 


-1/2 



The scalar doublets H u and H d are responsible for giving mass to the up and down quarks 
respectively. We refer to these doublets as the quark friendly doublets. Of the new doublets, 
the lepton friendly doublet Hi gives mass to the leptons, while the remaining inert doublet 
H does not couple to quarks or leptons. This Yukawa structure is enforced by a discrete 
Z2 symmetry, under which the superfields E, Hq, and Hi transform as X — > —X while all 
other fields remain unchanged. The most general superpotential respecting R-parity, gauge 
symmetry, and the Z2 symmetry is 

W = y u UQH u - y d DQH d - y e ELHi + Ji x H u H d + Ji 2 H Hi . (6) 

The Z2 symmetry is softly broken by the terms (n\H u Hi + $H§H d + h.c.) contained in 
the Higgs sector soft SUSY breaking potential Vs ft given by 
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The Higgs sector potential is given by the sum of the F-terms, D-terms, and Vsoft respectively 



V 



dW 



dHi 
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Soft- 
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Expanding the above expression results in 

V = m 2 u \H u \ 2 + m 2 d \H d \ 2 + m 2 \H \ 2 + m 2 \H e \ 2 + {n\H u H d + fi 2 2 H H e + ^H u H t + n\H Q H d + h.c.) 



+ 
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+ 



si 



|#d| 2 + |#o| 2 



\H, 



where m 2 = (l//^ 2 + M 2 ), m 2 d = (M 2 + fi 2 d ), m 2 = (|/i 2 | 2 + /j 2 ), m 2 = (\fi 2 \ 2 + /i 2 ), and a a 
(a = 1,2,3) are the Pauli matrices. To achieve spontaneous symmetry breaking, the Higgs 
doublets acquire the following vacuum expectation values (vevs): 

1 ( \ 1 / v d \ 1 ( 



v u J \ j i u„ 

We define t> 2 = v 2 + t>J + t>o +u| so that we have v 2 = AM z /(g 2 + gl) ~ (246 GeV) 2 . Between 
the quark friendly doublets we define the mixing angle tan/3 = v u /v d while between the 
lepton friendly and inert doublets we define the mixing angle tan fa = Vo/vg. We also define 
tana = v q /vi, where v 2 = v 2 + v d and v\ = v 2 + v 2 . These definitions allow us to express 
the individual vevs in terms of the Standard Model vev and the three mixing angles a, fa 
and fa 

v u = v sin a sin fa v d = v sin a cos fa v = v cos a sin/3^, vg — v cos a cos fa. (8) 

Each of the four complex Higgs doublets contains four real degrees of freedom, so there 
are a total of sixteen degrees of freedom. Three of these are eaten to give mass to the W ± 
and Z°, while those remaining result in a scalar mass spectrum that includes four neutral 
scalars, three pseudoscalars, and three charged pairs. From the scalar potential above, the 
mass matrices can be calculated. We parameterize them in terms of the gauge boson masses 
and the three mixing angles appearing in equation [HJ 

The neutral scalar mass matrix is M^ = 

/ M 2 -\M 2 s 2 a s w 

\M 2 s 2 a s w - n\ M 2 



\Mls 2a SfiSf }l -^M%S 2 aSpCi3 e ~ li\ \ 

-\Mls 2a cps Pt - nl \M\s 2a c fi Cf il 



\M 2 z s 2a spSp t -\M 2 z s 2a cps Pt - 14 Ml -\M 2 z c 2 a s We - \l\ 



\ -\M 2 z s 2a s fi c fil - /4 \M 2 z s 2a c p c^ -\M 2 z c 2 a s Wl - y? 2 



M 2 



where s x and c x are shorthand for sinx and cos a: respectively, and the diagonal terms are 
given by 

n\ cot/3 + n\ cot a (j? 
n\ tan p + fi\ cot a 
/i 2 cot fa + n\ tan a 



Ml 
Mf 



M| sin 2 a sin 2 + Ai, 
M| sin 2 a cos 2 + A 2 , 



Ml cos 2 a sin 2 /3 £ + A 3 , 



'2 ' 

r2 n „ P 2 



Ai 
A 2 

A3 



cos pi 
sin /3 

sin /?£ 
cos 8 



M\ cos 2 a cos 2 + A 4 , 
The pseudoscalar mass matrix is 



A 1 



/i 2 , tan /?£ + /i| tan a 



Mi 



( Ai /i 2 /x 2 \ 
Ai 2 A 2 /il 
//i A 3 a*1 



The charged scalar mass matrix is 



cos 8 
sin8 e 

sin P 
cos Pi 



(9) 



.U 2 . = M\ + AM 2 



(10) 



where 



AM 



^2/3 



1 c 2 

2 S a S 2P 



\ 



2 S 2a s P s p e 
\s2aSBCBt 



„2 2 
1 



C a C 2ft 



2 S 2a c P c Pi 



\s2 a SpSf3 e 
\s2aCpSp e 

+ S l C 2P 
l r 2 „ 



2 S 2aSpCfi i 

\s2aCpCi3 e 
l r 2 „ 



r 2 2 



S a C 2/3 / 



In Section 3.3 of |23j Gupta and Wells outline a procedure for obtaining an upper bound 
on the tree-level mass of the lightest neutral scalar, h, in the limit of large SUSY breaking 
masses (as compared to the Z-mass). The procedure consists of transforming the mass 
matrices into the so called "Runge basis," in which one doublet obtains all of the vev whi 



the others are orthogonal to one another. Details on the Runge basis can be found in 24]. 



In this basis all but one diagonal entry of the neutral scalar mass matrix grow large in the 
limit of large SUSY breaking masses. This entry acts as an upper bound on M\ since, for a 
positive definite matrix, the smallest eigenvalue is bounded above by the smallest diagonal 
entry. Their result holds in our case as well and results in the inequality 



Mh < M z I sin 2 a cos 2/3 + cos 2 a cos 2(3 e \ . 



(11) 
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Leading order radiative corrections to the Higgs masses will be important in constraining 
parameter space. As usual, the dominant contributions come from top quark loops, governed 
by the top quark Yukawa coupling. In this section we have written the neutral scalar mass 
matrix, Mjy, in the {u,d, 0,£} basis. Hence the 1-1 entry receives a correction from top 
quark loop diagrams given by 



where mi is the stop squark mass, which we take to be ~ 1 TeV. In addition to top quark 
loop corrections, other corrections are potentially significant because of the possibility of 
very large values for tan/3 and tan/3^. We therefore also consider the leading correction to 
the 2-2 and 4-4 entries of M^, which come from bottom quark loop diagrams and a tau loop 
diagram respectively. The 3-3 entry receives no correction since the inert doublet, H , does 
not couple to quarks or leptons. There are other sub-leading-log corrections to the masses, 
and these can contribute 5 — 10 GeV to the masses (see Ref. [25j for a detailed discussion). 

IV. CONSTRAINTS ON THE SUPERSYMMETRIC LEPTOPHILIC HIGGS 
MODEL 

In this section we outline the main constraints that limit the viable parameter space of 
the SLHM. The free parameters arising from the scalar sector consist of the four couplings 
fj%, A*3, and n\, which mix pairs of Higgs doublets in the scalar potential, as well as the 
three mixing angles tan a, tan /3, and tan/3^, which appear in equation [HI The constraints 
ari sing from the charged scalar sector are similar to those of the L2HDM, which is studied 
in [15|. Our interest therefore lies in the neutral sector. We find that LEP data and other 
constraints severely restrict the size of the allowable parameter space, but leave enough room 
to comfortably fit the model a lightest neutral scalar mass substantially less than 110 GeV. 

A. Yukawa Coupling Perturbativity 

The first constraints come from requiring that the Yukawa couplings remain perturbative. 
By demanding that each Yukawa coupling remains smaller than 4tt we obtain the following 




(12) 



S 



three inequalities 

13 2 , 

520 2 , (13) 
1235 2 . 

One can see that the top quark Yukawa coupling becomes non-perturbative for small values 
of tan a or tan B while the bottom quark Yukawa coupling does so for small values of tan a 
or large values of tan 0. In addition, the tau Yukawa coupling becomes non-perturbative for 
large values of tan a or tan Be. . 



tan 2 a J \ tan 2 8 J m 2 
1 + - — 5— 1 + tan 2 B < 



tan 2 a J \ J ml 
1 + tan 2 a 1 + tan 2 Be < 5- ~ 



m? 



B. Tree Level Unitarity 

Requiring perturbative unitarity of fermion anti-fermion scattering places upper bounds 
on the fermion masses. The unitarity condition that must be satisfied is |9R(a/)| < 1/2, 
where aj is the Jth partial wave amplitude in the partial wave expansion of the fermion 
anti-fermion scattering amplitude. The scattering we consider occurs by the exchange of a 
Higgs boson. We obtain bounds from imposing the unitarity condition on the J = partial 
wave amplitude, which is calculated from a sum over s- and t-channel helicity amplitudes in 



the high energy limit. The procedure is described in detail in 26], where contributions to the 
partial wave amplitudes are provided for a general model. These contributions depend on 
combinations of the vector and axial vector Yukawa couplings. For the SLHM the resultant 



bounds are found to be (see [26[] for a clear discussion) 



G F m\ 2 ■ 2 a 
=■ < sin asm B, 

4:77^2 



< sin ex cos B } (14) 



GfttiI 

G F m 2 2 2 a 
=■ < cos a COS Be. 

Here we have used the bounds obtained for third generation fermions as their larger masses 
yield the most stringent results. The unitarity constraint prevents very large values for 
tan B, capping it at around 300. Several combinations of tan a and tan /3 values on the order 
of several tenths are also eliminated. 
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C. The Anomalous Muon Magnetic Moment 



As in the Standard Model, the magnetic moment of the muon receives a contribution 
from the one-loop diagram formed by connecting the muon lines on a muon-muon-photon 
vertex with a neutral Higgs boson. Only the lightest neutral Higgs is relevant since the 
contribution goes as the square of the ratio between the muon and Higgs masses. For the 
SLHM the contribution is 

a K 2 K f 1 * 2 (2-- 
M 8ttV J z 2 + x 2 (1 

where x = M^/m^ and 



z) 



dz. 



(15) 



K 



\U A ' 2 



'41 



cos 2 a cos 2 /3f 

If the Higgs mass, Mh, is assumed to be the same in the SLHM and the Standard Model 
then the contribution to the muon's magnetic moment from a light scalar in the SLHM is 
simply its Standard Model value multiplied by A' 2 . The value of A 2 however, remains ^ 1 
across the entire spectrum of parameter space, even for very large values of tan a and tan /3e- 
A review on the anomalous muon ma gne tic moment is given by j^] while current results 
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291 ] . In our case the contribution is much too small 



and uncertainties can be found in 
to produce any bounds. 

In addition however, there is a two-loop Barr-Zee effect |30|, which is generally more 
significant than the one-loop contribution discussed above. The Barr-Zee effect occurs by 
connecting an internal Higgs to an internal photon through a massive fermion loop and is 
given by 3jJ, [32f] . We consider such effects with third generation fermions in the SLHM and 
find that the contribution to the muon magnetic moment is 

8U n f(xt) 2U 21 f(x b ) U 41 f(x T ) 



Aa, 



am U$i 



+ 



+ 



4:7T 3 v 2 cos ^ ^ 3 sin 2a sin (3 ' 3 sin 2a cos ' cos 2 a cos (3i 
where Xf — m 2 /M 2 and the function f(x) is given by 



(16) 



/(*) 



2z(l - z] 



In 



x 



dz. 



z(l — z) — x 

Though the contribution from the tau loop diagram is suppressed by m 2 /M|, it is enhanced 
for very large tan/3^. In following (33 1 we measure how well these contributions compare to 
experiment with the quantity 



Xa 



Afl SLHM 

6.8 x 10" 10 
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where 6.8 x 10~ 10 is the theoretical uncertainty for a M in the Standard Model (used because it 
is larger than the experimental uncertainty). The result is that, though larger than the one- 
loop contributions, the two-loop Barr-Zee effect contributions are still too small to provide 
significant constraints on the parameter space. 



D. LEP Higgs Search Data 

The largest source of constraints for the neutral sector of the SLHM consists of LEP's 
failure to discover a neutral Higgs boson. If the lightest neutral scalar's mass is too small, 
one would expect LEP to have seen it, whereas for a mass Mh > 114.4 GeV, LEP data 
becomes irrelevant and no bounds can be obtained 34|. The production mechanism at LEP 



is the Higgs-strahlung process e + e — > hZ, and thus if the coupling, gzzh, between the 
lightest neutral scalar and Z-pairs is sufficiently small, the scalar's non-discovery at LEP 



can be explained 
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In addition, there is an effect which suppresses the sensitivity with which the experimental 
results may be applied to constrain models beyond the Standard Model 
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Bounds 



from LEP were produced under the assumption that the Higgs boson decays exclusively 
into bb pairs or exclusively into r + r~ pairs. LEP has provided a bound on the quantity 
BR(/i — > XX)£ 2 for X = b and X = r, where £ is the ratio of the ZZh coupling in a model 
to that of the Standard model i.e. £ = gzzh/ 1 9^ih- We find the value of £ 2 in the SLHM to 
be 



e 



U\\ sin a sin (3 + U 2 i sin a cos (3 + U S i cos a sin /3e + f/41 cos a cos j3g 



(17) 



We will employ both of these bounds to exclude regions of parameter space in the SLHM. 
Naively, one expects BR(/i — > bb) to approach unity when tan/3 is large and tana, tan/^ 
are small since in that case the down-type quark Yukawa couplings are doubly enhanced 
while the lepton Yukawa couplings remains small. On the other hand, when tana and 
tan/3^ are large while tan/3 is small, the lepton Yukawa couplings are enhanced and the 
down-type quark Yukawa couplings remain small, resulting in an increase in the branching 
ratio BR(/i r+r~). 

Since in the interesting region of parameter space, the ZZh and WWh couplings are 
small, we can approximate the total decay width as simply T(h — > bb) + T(h — > t + t~). 
The two branching ratios for the SLHM can therefore be conveniently expressed as BR(/i — > 
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bb) = 1/(1 + k) and BR(A t + t~) = k/(1 + k), where k = T(h T + T~)/T(h -»■ 66). The 
variable k is straightforward to calculate and is given by 

4m 2x3/2 



3m 2 



2 cos /3 
tan a- 



cos 2 f3e 



U. 



41 



21 



M, 2 



M 2 - 4m 2 



(18) 



where the Uij are entries of the 4x4 diagonalizing matrix defined by WM^U = M|j . 
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FIG. 1: The colored regions illustrate the allowed points in the tana, tan/3, tan/3^ parameter 
space. Each region is a slice of constant tan fig in the tana x tan/3 plane. The values of fii, fi2, 
/^3, and are fixed at 200, 250, 300, and 100 GeV respectively, but changing fj,i and/or ^3 has 
relatively little effect. Increasing [12 and/or /U4 shrinks the above space. Increasing tan/3^ enlarges 
the size of the allowed space quite rapidly until around tan/3^ « 8, when the space stops enlarging 
and begins to slowly shrink - this can be seen in Figure 



We have numerically scanned through parameter space, calculating the values of BR(/i — » 
&F)£ 2 ,BR(/i — > t + t~)£ 2 , and in the SLHM. Those points in parameter space for which 
either BR(/i — > bb)^ 2 or BR(/i — > t + t~)£ 2 is greater than its LEP bound at the corresponding 
value of Mh are excluded. By imposing these two LEP bounds as well as the perturbativity 
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requirements of Section IIV Al and the unitarity requirements of Section IIVB| we are able to 
exclude substantial regions of the model's parameter space. In Figures [1] and [2] the allowed 
region of the three-dimensional parameter space for the variables tana, tan/3, and tan/3^ 
is shown. For these plots the values of fii, /12, (13, and /14 have been fixed at 200,250,300, 
and 100 GeV respectively. The plots depict several sections of viable parameter space in the 
tana x tan/3 plane, each being a slice of constant tan/3^. As tan/3^ varies over its allowed 
range, one can see how the sections grow in area, change shape, and eventually shrink back 
away. 
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FIG. 2: A continuation of figured] for larger values of tan/?£. As tan/3^ increases beyond 80, the 
space very slowly shrinks into an extremely thin sliver of possible tana values centered near 2; it 
finally disappears completely at tan « 350. 

Though the values of fii, /i 2 , /x 3 , and /x 4 are fixed, the size and shape of the allowed 
parameter space remains largely unchanged when fii and /X3 are allowed to vary between 
50 and 1000 GeV. Their values are consequentially relatively unconstrained. Increasing the 
value of /i4 however, has the effect of sharply cutting down on the size of the allowed region 
of parameter space. So too does increasing /z 2 , though to a slightly lesser degree. Merely 
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increasing /X4 to 200 GeV results in a drastically smaller allowed region than that shown 
in Figure [T]and completely eliminates the regions corresponding to tan (3g values of 5.3 and 
5.5. The other regions are compressed so that 3 < tana < 20 and 50 < tan/3 < 290, while 
their overall shape remains the same. Enlarging either /i 2 or /i 4 further rapidly shrinks the 
allowed space away until it vanishes completely. 

Figure [3] plots an assortment of possible BR(/i — > bb)^ 2 values as a function of the lightest 
neutral scalar mass Mh- Each value plotted corresponds to some point in the allowed region 
of parameter space. The LEP curve is shown in blue. For very large values of tan/3^, the 
curves continue down to approximately 25 GeV, with the value of BR(/i — > bb)l; 2 becoming 
extremely small. We see that Higgs bosons below 114.4 GeV are certainly allowed, but 
below approximately 90 GeV their couplings to vector bosons become negligible, making 
detection through vector boson fusion or Higgs-strahlung off a vector boson impossible. The 
analogous result for BR(/i — > t + t~) is plotted in Figure HJ with similar conclusions. 
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FIG. 3: Various values of the quantity BR(/i — > bb)^ 2 plotted as a function of the lightest neutral 
scalar mass Mh- The plotted values correspond to a uniform sampling of points within the allowed 
regions of the tan a x tan /3 plane for the different values of tan j3 e u that are plotted in figures [T] 



and The LEP bound of reference 



361 ] is shown in blue. 
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neutral scalar mass M^- 

V. PHENOMENOLOGY 

In this section we discuss the possibility of detecting a supersymmetric leptophilic Higgs. 
We have focused on the neutral sector, as the charged sector strongly resembles the non- 



SUSY leptophilic scenario covered in [15] . The quantity of importance to the decay of the 
lightest neutral scalar is the ratio k = BR(/i — > t + t~ )/BR(/i — > bb), which is given by 
equation [18] in Section IIVDI 

For the region of parameter space discussed in the previous section, we have shown various 
values of k in Figure [5j For Higgs bosons near 114.4 GeV, the allowed value of k approaches 
its Standard Model value of approximately 0.1. However, for lighter Higgs bosons, k is much 
bigger, approaching unity for Higgs masses below 100 GeV. 

We see that in this model, the Higgs can be relatively light, and will have a much larger 
branching ratio to r + r~ than in the Standard Model. In order to detect the Higgs at the 
Tevatron or the LHC, however, one also must consider the production rate. As we have seen, 
for Higgs bosons below 90 GeV, the ZZh and WWh couplings are quite small, and thus 
Higgs-strahlung is negligible. What about gluon fusion, which is the primary production 
mechanism for a light Higgs? Here, one must include both top and bottom loops, and the 
coupling to the Higgs will be different. We find that the ratio of the gluon fusion cross 
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section to that of the Standard Model is 



Q" SLHM 



u 
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A(m b ) U, 
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sin a sin /3 A(m t ) sin a cos /3 



(19) 



and this is plotted in Figure M for various parameters. The function A(m,f) is given by 
A(m,f) = 2[xj + (xf — xj 2 , where Xf = M^/Am 2 and f(x/) is given by equation 

2.47 in [40]. For much of parameter space, the gluon fusion rate is also very small, making 
Higgs detection extremely difficult. In the Standard Model, the only other production 
mechanism that doesn't involve gluon fusion or the WWh or ZZh vertex is Higgs-strahlung 
off a top quark. That is difficult in the Standard Model, and in this model is even weaker 
since the top quark Yukawa coupling is smaller. One can think about Higgs-strahlung off a 
tau, but this is likely to be swamped by backgrounds. 

In any event, this is just a specific model. One might have other possibilities for Higgs 
production, such as production in the decay of one of the charged Higgs bosons in the model, 
or production through supersymmetric particles. In both of these scenarios, the production 
rate would depend on many additional parameters. Thus, experimenters should look for 
Higgs bosons in the 75 — 110 GeV range with a substantially enhanced coupling to r pairs 
(below 75 GeV, a very small sliver of parameter space does remain). A study of r pair 
detection in leptophilic Higgs decays at the LHC was carried out in Ref. 4l). Since they 
did not consider the supersymmetric version, they concentrated on Higgs in the 100 — 160 
GeV mass range, and gluon fusion production was not particularly suppressed, as it is here. 
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FIG. 6: Logplot of the ratio of the production cross section of the lightest neutral scalar by gluon 
fusion in the SLHM to the Standard Model. 

They also focussed on models with dark matter candidates (usually involving an additional 
singlet or an additional inert doublet). Nonetheless, their techniques show that detection of 
a Higgs decay into r pairs is feasible in the early stages at the LHC. At the Tevatron, CDF 
and DO did explicitly search for Higgs decays to r pairs 42j, but did not consider Higgs 
masses below 90 GeV 

Throughout this analysis, we have ignored the effects of the heavier neutral Higgs scalars. 
Consider the second lightest neutral scalar, rj. As we scan the entire allowed parameter 
space, we find that the r] always appears to be very close to 110 GeV. This may not be too 
surprising. Imagine that there was no mixing at all between the quarkophilic and leptophilic 
Higgs sectors. Then each sector would have a similar mass matrix to that of the MSSM 
(although with smaller overall vevs), and thus one would find two relatively light Higgs. 
Mixing can't be eliminated, of course, due to D-terms, but it is not surprising that there 
are two relatively light scalars in the model. In the region of parameter space in which 
the couplings of the h to the gauge bosons is severely suppressed, however, the couplings 
of the rj will not be, and thus the rj will be similar to the Standard Model Higgs. Given 
the uncertainty in our calculations, including the effects of non-leading-log and higher order 
corrections to the masses, it is premature to conclude that the current LEP bounds would 
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rule out this 110 GeV Higgs, but an increase of just a few GeV in the current lower bound 
on the Standard Model Higgs would rule out this model. 

In the region of parameter space of interest, the h and i] are primarily linear combinations 
of H and H u , with small admixtures of H& and He. Nonetheless, the ratios of vacuum 
expectation values are large enough that the dominant decay of the h, for example, is 
primarily into r's and 6's through these small admixtures. The two heaviest Higgs bosons 
are each almost entirely Hd and He , respectively, with little mixing. 

Consider these two heavier Higgs bosons, Hi and H 2 . Since the coupling of the 77, in 
the region of interest, to Z-pairs is very close to that of the Standard Model, then the 
fact that the sum of the squares of the Higgs couplings to Z-pairs must equal the square 
of the Standard Model coupling implies that the coupling of Hi and H 2 with W, Z-pairs 
is negligible. We have confirmed this numerically. Another way to say this is that the 
narrow window of parameter space forces the direction of the vacuum expectation value to 
be almost entirely in the 77 direction, leaving little room for vev-dependent couplings of the 
other neutral Higgs. This will also cause a suppression in the Hihh and H 2 hh couplings. The 
Hi and H 2 will thus be both Higgs-phobic and gauge-phobic and will only decay into fermion 
pairs. One of the two, Hi, will decay almost entirely into bb, and the other, H 2 , will decay 
almost entirely into t + t~. This leads to interesting phenomenological consequences. The 
Hi can be copiously produced through gluon fusion (through its coupling to the 6-quark), 
and its dominant decay into bb will be quite dramatic. The H 2 would be a heavy Higgs 
boson that decays entirely into r pairs. However, gluon fusion occurs at a small rate, and 
thus production through heavier particles or supersymmetric partners would be necessary. 
This possibility is currently under investigation. 

VI. CONCLUSION 

In this work, we have studied the Higgs sector of the supersymmetric version of lep- 
tophilic models. The model contains four Higgs doublets, which couple to the up quarks, 
down quarks, charged leptons and no fermions, respectively. The Higgs sector, as in all 
supersymmetric models, is tightly constrained. We consider constraints from perturbativ- 
ity, unitarity, the muon anomalous magnetic moment and we also impose constraints from 
experimental searches at LEP. 
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We find that in most of parameter space, the lightest Higgs, h, has a mass between 75 
and 110 GeV (with a very small sliver of parameter space giving smaller masses). For lighter 
values of the mass, the decay branching ratio into r pairs is substantial, and can even be 
the dominant decay mode. This would lead to some spectacular signatures at the Tevatron 
and the LHC. However, the conventional production mechanisms, such as W-fusion, Higgs- 
strahlung and gluon fusion are suppressed in this region of parameter space. 

The second lightest Higgs, rj, has a mass throughout the allowed parameter space of 
approximately 110 GeV. Its production cross section is not as strongly suppressed, and would 
appear similar to a Standard Model Higgs. The remaining two neutral scalars are typically 
heavier, are gauge-phobic and Higgs-phobic, and would decay into fermions. One decays 
almost entirely into bb and would be copiously produced through gluon fusion. The other 
decays almost entirely into r + r~, but conventional production mechanisms are suppressed. 

There are also three charged scalars and three pseudoscalars in the model. We do not ex- 
pect the phenomenology to differ substantially from the detailed analysis of Logan and 
MacLennan[15], who used MSSM parameters to constrain their parameter space (even 
though the model was not supersymmetric), and thus there would only be 0(1) changes 
in their results due to mixing angles. Exploration of the supersymmetric particles in the 
model are currently under investigation. 

We thank Heather Logan and Reinard Primulando for useful discussions. This work was 
supported by the National Science Foundation PHY-0755262. 
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